Unconventional Fermi surface instabilities in the Kagome Hubbard Model 
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We investigate the competing Fermi surface instabilities in the Kagome tight-binding model. 
Specifically, we consider onsite and short-range Hubbard interactions in the vicinity of van Hove 
filling of the dispersive Kagome bands where the Fermiology promotes the joint effect of enlarged 
density of states and nesting. The sublattice interference mechanism [Kiesel and Thomale, Phys. 
Rev. B Rapid Comm., in press.] allows us to explain the intricate interplay between ferromagnetic 
fluctuations and other ordering tendencies. On the basis of functional renormalization group used 
to obtain an adequate low-energy theory description, we discover finite angular momentum spin 
and charge density wave order, a two-fold degenerate d-wave Pomeranchuk instability, and /-wave 
superconductivity away from van Hove filling. Together, this makes the Kagome Hubbard model 
the prototypical scenario for several unconventional Fermi surface instabilities. 



PACS numbers: 71.10.Fd,71.10.-w,64.60.ae,74.20.Mn 

Introduction. The interplay of Fermiology and inter- 
actions gives rise to a plethora of ordering phenomena 
in two-dimensional electron systems. Starting from an 
itinerant electron picture, the density of states (DOS) at 
low energies around the Fermi level as well as nesting fea- 
tures of the Fermi surface are the relevant parameters of 
the kinetic theory. In the limit of weak interactions im- 
posed on the non-interacting electrons where a perturba- 
tive treatment is asymptotically exact, superconducting 
order as a phase-coherent superposition of Cooper pairs is 
the generically encountered Fermi surface instability pQ. 
This situation can be changed in various different ways: 
By enhancing the interaction scale or via nesting and en- 
larged DOS at the Fermi level, order due to condensation 
of particle-hole pairs can become competitive and even 
favorable to superconductivity. Prominent examples in- 
clude magnetic (charge) order via a spin (charge) density 
wave which can induce superconductivity as a function 
of doping or pressure. 

Since the discovery of the cuprates, it has widely been 
appreciated that electronically mediated interactions in 
particular, can give preference to electron condensates 
with finite angular momentum of the condensing pairs. 
It implies that the pairs which condense form at finite dis- 
tance, so as to minimize Coulomb repulsion in the case 
of a Cooper pair. This yields a momentum dependence 
of the associated mean field order parameter which can 
impose nodes on the Fermi surface in the ordered phase. 
In the case of a particle-hole condensate with opposite 
charge of the pair constituents, it is less generic that 
a pair of finite angular momentum should be energeti- 
cally preferable. As another difference to particle-particle 
pairs, the orbital angular momentum of a particle-hole 
pair does not unambiguously determine the spin of the 
pair to be singlet or triplet, which implies an even richer 
variety of possible orderings [2H1] . While this is short of a 
complete characterization, one important factor in favor- 
ing such orders is given by a tuned arrangement of longer 



range interactions [3J, a direction which has become re- 
cently accessible experimentally with sufficient tunability 
in dipolar fermion models [5] . Still, until today the main 
challenge in theory has been to find bare models of inter- 
acting electrons where these phases can be found as the 
natural ordered state at low energies. 

In this Letter, we propose and analyze the Kagome 
Hubbard model (KHM) as a prototypical microscopic 
model to realize certain kinds of such unconventional 
Fermi surface instabilities. The Kagome lattice [B] pos- 
sesses a minimal three-band model due to three sites per 
unit cell (Fig. [iji). For the KHM, the three sublattices 
have a fundamental impact on the preferred electronic 
many-body state at all regimes of coupling strength. In 
the strong coupling limit at half filling, the Kagome spin 
model exhibits strong quantum spin fluctuations and has 
become a primary candidate for quantum spin disordered 
phases [THH]- At intermediate coupling, comparatively 
little is known about the electronic quantum phases. At 
fillings up to the flat band (Fig. [T^,), ferromagnetism has 
been proposed on the basis of Stoner's criterion |10) . 

For the dispersive band fillings which we address in this 
work, the scenario is more complicated and involves the 
interplay of nesting, Fermi level DOS, and interactions. 
This is likewise indicated from recent studies at infinites- 
imal coupling where the sublattice interference mecha- 
nism has been developed as a key property to understand 
the KHM [11] . Out of the search for microscopic mate- 
rial scenarios of a Kagome lattice model at intermediate 
coupling, the Herbertsmithites such as ZnCu3(OH) 6 Cl2 
appear as an important class of candidates, while it is 
still hard to judge from neutron experiments how rele- 
vant charge fluctuations are and how the system behaves 
upon doping [12]. Promising alternative routes start to 
emerge in optical Kagome lattices of ultra-cold fermionic 
atomic gases such as for the isotopes 6 Li and 40 K [T3] . 

Main results. As a function of local and nearest neigh- 
bor Hubbard coupling Uq and Ui, we find a rich phase 
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FIG. 1. (Color online). Fermiology of the Kagome tight- 
binding model, (a) Band structure where the shaded region 
is the vicinity around van Hove filling (VHF) n = 5/12 with 
an enlarged density of states (inset), (b) Fermi surface at 
VHF. An N — 96 patch discretization of the Brillouin zone 
is given along with the change of dominant sublattice occu- 
pation (red, blue, and green dots) on the Fermi surface. The 
nesting vectors Qi,2,3 connect different sublattice states at 
the Fermi level. The grey dashed line hints the Fermi surface 
at n = 5/12 + 0.02. 



diagram of the Kagome Hubbard model (KHM) at and 
around van Hove filling (VHF) which is summarized in 
Fig. [2] Right at VHF where the Fermi level DOS is 
maximally enhanced and nesting features of the Fermi 
surface are strongest (Fig. [Ij}) , the system promotes fer- 
romagnetism (FM) for dominant Uq. This is a conse- 
quence of the suppression of local Hubbard matrix ele- 
ments due to sublattice interference which otherwise 
would give rise to spin density wave or superconducting 
order. As U\ is enhanced, we discover a p-wave (L = 1) 
charge bond order (cBO) and spin bond order (sBO) 
phase (Fig. [3]). This is again motivated by the sublat- 
tice interference which suppresses the energy gain of a 
zero angular momentum particle-hole condensate. For 
dominant U±, we find a d-wave (L = 2) Pomeranchuk 
instability (PI) which is two-fold degenerate due to the 
associated two-dimensional irreducible point group rep- 
resentation of the Kagome lattice, i.e. the E% element of 
the Cq v symmetry group. This explains why the nematic 
phase resulting from there can break the rotation symme- 
try of the Kagome lattice in different ways, which leads 
to different distortions of the Fermi surface (Fig. [4]) . As 
we deviate from VHF, nesting effects get reduced. While 
the PI (FM) phase still persist for dominant U\ (Uq), 
/-wave superconductivity (/-SC) emerges in the inter- 
mediate regime as a consequence of longer range interac- 
tions and ferromagnetic fluctuations which promote spin 
alignment of the Cooper pairs (Fig. [5]) . 

Kagome Hubbard model. We consider the Hamiltonian 

h = Yl E ( c ls- + h - c -) + /* E n ^ 

+Uq ^ »i,t"i,4. + "y n>i,anj,<r' , (1) 

where Uq denotes the local and U\ the nearest neighbor 
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FIG. 2. (Color online). Phase diagram of the Uq-U\ Kagome 
Hubbard model at (a) VHF n = 5/12 and (b) away from VHF 
at n — 5/12 + 0.02. The shaded areas for small Uq and U\ 
indicate regimes where critical ordering scales were too small 
to be determined. In (a) ferromagnetism (FM) is found for 
large Uq, along with a d-wave Pomeranchuk instability (PI) 
for large Ui as well as intermediate spin bond order (sBO) and 
charge bond order (cBO) phases. In (b) the PI and FM phase 
persist away from VHF, while the density wave orders disap- 
pear in favor of an /-wave superconducting (/-SC) phase. 



Hubbard term, /i is the chemical potential, and the tight 
binding hopping term has been set to unity. In what fol- 
lows, we adjust the chemical potential such that we are 
around VHF n = 5/12, and set the Fermi level at VHF 
to zero energy (Fig. [T^,). There, the Fermi surface pos- 
sesses a hexagonal form where the van Hove points are 
located at the M point in the Brillouin zone (Fig. [IJd) . 
It follows that there are three important nesting vectors 
Qi = tt(-|,-^), Q 2 = tt(1,0), and Q 3 = tt(-§, ^), 
with the lattice constant between two adjacent sites set 
to unity. Away from VHF, the Fermi surface is rounded 
(gray dashed line in Fig.[ljD) and nesting is reduced. The 
nesting vectors Qi,2,3 connect parts of the Fermi sur- 
face whose states are dominated by different sublattices 
(Fig. [IJd) . As a consequence, by transforming the local 
Hubbard term Uq into band space, the matrix elements 
along the nesting vectors are suppressed which we call 
sublattice interference [TTj. If this interference mech- 
anism were absent or neglected, the short-range KHM 
would take on phases dictated only by the nesting vectors 
Qi,2,3, such as a spin density wave state, and d-wave su- 
perconductivity below some finite coupling strength [14j . 
Instead, the sublattice interference enhances the rele- 
vance of ferromagnetic fluctuations stemming from the 
local Hubbard term, and also promotes the relevance of 
the nearest neighbor Hubbard coupling. 

N -patch functional renormalization group. We employ 
the functional renormalization group (FRG) to obtain an 
effective low energy description of the bare model in ([!]) , 
which has proved suitable in a variety of interacting two- 
dimensional electron systems [TS]. We study how the 2- 
particle vertex evolves under integrating out high-energy 
fermionic modes along a temperature-flow cutoff scheme: 

V A (fci,fe2;fc3,fc4)4 4 , s 4 3 , s - c fc 2 ,s c fe 1 ,,' where A is the IR 
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FIG. 3. (Color online). Real space patterns of L = 1 density 
wave orders with (a) charge bond order (cBO, S = 0) and (b) 
spin bond order (sBO, S — 1). Both orders are superpositions 
of the the one-dimensional alternating bond orders induced 
along the momentum transfer Qi,2,3, yielding a 12-site unit 
cell. Enhanced charge density bonds in (a) are indicated by 
thick black lines, f (4) bonds in (b) by thick blue (red) lines. 



cutoff approaching the Fermi surface within the flow, fci, 
hi (fea, ki) denote the ingoing (outgoing) fermionic mo- 
menta, and s, s take on opposite spin orientations. This 
is sufficient because for a spin-rotation invariant model, 
the S z = sector of the scattering vertex allows us 
to extract the (triplet) singlet channel by (anti-)sym- 
metrization of the vertex. We neglect the self energy 
corrections along the flow and discretize the k's to rep- 
resent specific patches in the Brillouin zone. In Fig. [TJd, 
we have depicted such a patching scheme for TV = 96 
patches. The phase diagrams in Fig. [2] have been ob- 
tained with this level of discretization. In order to assure 
that the discretization yields converged results, we have 
employed supercomputer facilities to compute selected 
points of the phase diagram for up to N — 384 patches, 
which corresponds to solving a 5.7 x 10 7 -dimensional sys- 
tem of integro-differential equations. The leading diverg- 
ing channel of the 2-particle vertex signals the occurrence 
of a Fermi surface instability. It is then decomposed into 
eigenmodes to obtain the associated form factor [16] . 

Ferromagnetism (FM). As alluded to above, the sub- 
lattice mechanism suppresses strong finite momentum 
scattering channels originating from Uq. When the local 
Hubbard interaction is dominant, the large Fermi level 
DOS drives ferromagnetic fluctuations without any other 
competing channel, which thus explains the formation of 
FM order. This holds at and around VHF (Fig. [2). The 
order parameter reads O fm = J2k,i,s.s'( c kis a '^'C kls ,), 
where er denotes the vector of Pauli matrices. In addi- 
tion, the propensity towards spin alignment from the fer- 
romagnetic fluctuation background at high energies pro- 
vides further bias for the sBO phase at VHF and the 
/-wave SC phase away from VHF. 

Spin and charge bond order (sBO and cBO). As U\ 
is enhanced at VHF, the nesting vectors Qj.,2,3 become 
important again. This is because the sublattice interfer- 
ence does not apply to nearest neighbor Hubbard inter- 
actions which connect different Kagome sublattices. As 



a consequence, the KHM exhibits cBO also known as the 
Peierls phase which is characterized as a L = 1, 5 = 
particle-hole condensate as well as its spinful condensate 
counterpart L = 1, S = 1, i.e., the sBO phase (Fig. |3j. 
The cBO shows an alternating sequence of bonds with 
enhanced and reduced charge. Accordingly, the sBO ex- 
hibits alternating bonds of enhanced j" spin occupancy 
and enhanced J. spin occupancy. Both cBO and sBO 
can be deconstructed into separate bond orders along the 
three individual directions in the Kagome lattice, which, 
recast in momentum space, are exactly Qi,2,3- For sBO, 
the order parameter is given by 

h,s,s f ^ ' 

l,m.n 

(2) 

where £i mn is the Levi-Civita tensor. Eq.[2]takes the sim- 
ilar form for cBO, where cr ss i is replaced by the identity 
matrix l ss >. There are in principle three independent 
bond order mean fields associated with these directions. 
As they are independent and degenerate, however, the 
system gains energy by forming all of them at the same 
time, which results in a 12-site unit depicted in Fig. [3ji 
for cBO and in Fig.[3}D for sBO. We checked in a 12-band 
mean field free energy analysis that the system linearly 
gains energy from forming the individual mean fields, i.e. 
the ordering along the individual bond directions is in- 
dependent. 

Pomeranchuk instability (PI). The PI forms for dom- 
inant U\ at and around VHF. The singularity in the 
particle-hole channel is located at momentum transfer 
Q = 0, implying that it does not break translational 
symmetries, but instead drives the system into a nematic 
phase as it breaks the lattice rotational symmetries of the 
KHM. The particle-hole condensate is characterized by 
L = 2, S = 0. According to the irreducible lattice rep- 
resentations of the Kagome lattice, a ci-wave instability 
necessitates a two-fold degenerate d x 2_ y 2 , d xy subspace of 
solutions [T7] . These form factors are depicted in Fig. [4^,. 
It turns out that the most dominant harmonics in this 
symmetry sector relate to the third nearest neighbor, i.e. 
equal sublattice particle-hole pairing. It is revealing to in- 
vestigate what kind of Fermi surface distortion can result 
from such an instability. As stated before, the PI effec- 
tively generates a third nearest neighbor hybridization. 
Its lattice rotational character, however, is not uniquely 
specified because we could in principle form any real su- 
perposition of the d x 2_ y 2 1 d X y solutions |18j . i.e. 

fd ;c 2_ y 2 ( k ) = cos ( 2k x) ~ cos (k x ) cos (yZky^j , 

fd xy (fe) = VSsin (k x ) sin (Vsky) . (3) 

Within our approach, we are short of an answer which 
linear combination is energetically preferable, since the 
lack of self energy damping results in an unbound energy 
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FIG. 4. (Color online). Pomeranchuk instability (PI) in the regime U\ > Uo. The Q = particle-hole pairing does not break 
translational symmetries and is dominated by third nearest neighbor pairing of the same Kagome sublattice which follows from 
the form factors of Eq. [3] shown in (a). The PI lies in the L = 2 sector (d-wave type) and is two-fold degenerate. It yields a 
distortion of the Fermi surface which can break the rotational symmetry in different ways such as to Ci v in (b) or C2 in (c). 



gain from Fermi surface distortion. This ambiguity is 
also visible in the order parameter 
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(k) + Bf dx (k)), (4) 



where A and B specify the superposition. We note, how- 
ever, that the remainder rotational group down to which 
the PI can establish a nematic phase depends on the cho- 
sen linear combination. In Fig. [4Jd, the remainder group 
is C% v and the Fermi surface is shifted away from all van 
Hove points, while for example in Fig. Hb one van Hove 
point remains unaltered and the remainder group is GV 
j '-wave superconductivity (f-SC). Away from VHF for 
intermediate Uq and TJ\ we find an /-wave, i.e. L = 3, 
5 = 1 Cooper pair condensate in the KHM. Fig. [5^, shows 
the form factor associated with the order and Fig. [5)3 
depicts the real space amplitude pattern which is essen- 
tial to identify the dominant harmonics in the /-wave 
symmetry sector, i.e. the typical pairing distance |17j . 
We find that SC pairing forms between second near- 
est neighbors to avoid the short-range Hubbard inter- 
action. Denoting the three sublattices by a, /3, and 7, 
we can decompose the pairing into three different form 
factor contributions given by f a ,p(k) = sin(^k x + ^k y ), 

//3, 7 (fe) = sind*^ - Y-hy), and /« )7 (fc) = sin(v / 3fe y ), 
where f m , n = fn,m and f m<m = 0. Together, they form 
the order parameter 
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fm.n {k) , (5) 



which can be rotated within the triplet state to the S z 
■ 1 sectors. Such a multi-sublattice scenario is similar to 
a multi-orbital superconductor where the form factor can 
be similarly decomposed into different orbital pairs [T§] . 

Note added. Upon completing the manuscript, we be- 
came aware of an independent work that investigates 




FIG. 5. (Color online), /-wave superconductivity (/-SC) 
away from VHF for Uo ~ U\. (a) /-wave pairing form factor 
(black solid line) between the a and j3 sublattice. The signal 
should be small for the 7-dominated domains on the Fermi 
surface (red patches 1 — 16 and 48 — 64) . The a and /? regimes 
of the form factor are best fit via second nearest neighbor 
pairing harmonics (dotted lines). The real space amplitudes 
of the a-/3 pairing is shown in (b). The total /-wave signal is 
built by adding the analogous a-7 and f3-j components. 



the Kagome Hubbard model strictly at van Hove filling 
through singular mode FRG, which is a complementary 
approach to iV-patch FRG with enhanced radial resolu- 
tion of the RG flow [2D]. The discrepancies in terms of 
qualitative results could trace back to assumptions made 
in the singular mode FRG approach where the vertices 
are reduced to their bosonic transverse momentum de- 
pendence in the respective channel |21j . This approxima- 
tion might be less valid when ordering of finite distance 
pairs emerges yielding important dependencies of all ver- 
tex momenta. This is the generic case for most Kagome 
lattice orders we find such as PI (third nearest neighbor) 
or /-wave SC (second nearest neighbor). 
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